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We consider the chiral limit of QCD subjected to an imaginary isospin chemical potential. In the 
e-regime of the theory we can perform precise analytical calculations based on the zero-momentum 
Goldstone modes in the low-energy effective theory. We present results for the spectral correlation 
functions of the associated Dirac operators. 
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I. INTRODUCTION 

An isospin chemical potential provides a way to "twist" 
the usual Dirac operator in two different directions. A 
real isospin chemical potential [l| gives a fermion deter- 
minant that is real and positive, and thus amenable to nu- 
merical simulations; it does not, however, preserve anti- 
hermiticity of the Dirac operator itself. An imaginary 
value of the isospin chemical potential, on the other hand, 
gives massless Dirac operators that are anti-hermitian; 
their eigenvalues thus lie on the imaginary axis instead 
of spreading out into the complex plane. This makes 
imaginary isospin chemical potential a useful parameter 
for deformation of the Dirac eigenvalue spectrum. 

Recently we have noted 0,|3 that a particular spectral 
two-point correlation function of Dirac operator eigen- 
values near the origin yields a direct way of determining 
the pion decay constant, F-^, from lattice gauge theory 
simulations. An alternative proposal, also using imagi- 
nary isospin chemical potential, is to use the distortion 
of the mass-dependent chiral condensate^ [j. For the 
quenched theory it has also been shown that one can 
use an ordinary (real) baryon chemical potential to ex- 
tract the pion decay constant from lattice distributions 
of the Dirac eigenvalues 0,13 ■ For gauge group SU(3) 
and quarks in the fundamental representation, the Dirac 
operator spectrum with real baryon chemical potential is 
complex, and the theory with dynamical quarks is dif- 



ficult to simulate directly due to a complex fermion de- 
terminant. For this reason imaginary isospin chemical 
potential is a more convenient strategy. 

For a given non-Abelian gauge potential Afj^{x) we 
study the two Dirac operators 



i?+vi"^ ^ m^) + *wso7o]^i"^ = *aV"VA"^ (1) 



and 



^ The twisted boundary conditions used there are equivalent to an 
imaginary isospin chemical potential. 



i?_^L"> ^ [Ip{A) - zMiso7o]V'i"^ = *AL"Vi"\ (2) 
where /iigo is real. Both operators D± are anti-hermitian, 

(n) 

and the eigenvalues Aj_ therefore lie on the real line. 
An imaginary isospin chemical potential can be viewed 
as an external constant Abelian gauge potential ^o that 
couples to the u and d quark with opposite charges. 

Much work has gone into understandin g gaug e theo- 
ries at real isospin chemical potential 0, H S U^ IHl IHi 
Il3 . Il4| . often in terms of the effective low-energy the- 
ory. Here we consider the effective theory, a chiral La- 
grangian, in the presence of an imaginary isospin chemi- 
cal potential. With the usual pattern of spontaneous chi- 
ral symmetry breaking for two light flavors, the theory is 
described by a Lie-group valued field U{x) € SU(2). 

Our focus will be on the so-called e-regime of QCD, 
where the chiral Lagrangian is treated as a perturbative 
expansion around the zero-momentum modes in a finite 
volume V |15J| . Roughly speaking, we are dealing with 
an expansion in 1/ L (where L is the linear extent of the 
finite volume) rather than the usual expansion in a small 
momentum p. There is a well known and intriguing con- 
nection between this regime and a universal limit of Ran- 
dom Matrix Theory [li [Hill [H ^ HI HI , but here 



we will stay entirely within the framework of the effective 
chiral Lagrangians. 

When the chemical potential is included, the power 
counting of the e-expansion must be reconsidered. One 
factorizes the field as U{x) — Uoexp[i\/2(j){x)/F.^], where 
Uq is the zero momentum part that will be treated ex- 
actly, and (f>{x) represents the fluctuation fields (without 
zero modes). It turns out that the naive guess provides 
a consistent counting: To leading order one keeps only 
the static modes in the path integral, while the fluctua- 
tion degrees of freedom decouple. This is not completely 
obvious at first glance, but can be seen as follows. The 
coupling to the imaginary isospin chemical potential in 
the chiral Lagrangian is dictated by the way it couples 
at the quark level [Eqs. ^ and |(2Jl]. Vector sources of 
that kind give rise to a covariant time derivative in the 
effective SU(2) Lagrangian dH^, 

doU{x) ^ VoU{x)^^oU{x)-^^iiso[<y3,U{x)l (3) 

where a^ is the usual Pauli matrix. The leading-order 
terms in the effective Lagrangian then read 

F^ 



C = ^Tt [VoU{x)VQU\x) + ^^U{x)^iU''{x)] 
-|Tr [MU^x)+M'<U{x)] , 



(4) 



S is the chiral condensate. When we expand 

U{x) = C/o [l + tV2(j>{x)/F^ + •••], (5) 

this produces the usual kinetic term for (t>{x). Let us 
recall the power counting in the e-expansion |l5| : We 
assume TOtt ^ p^ — 0{e^) while (f){x) ^ 1/L = 0{e), and 
a consistent power counting for the /iiso-term is /Ziso = 
O(e^). Indeed, when we expand the covariant derivative 
(O using Eq. (0 the leading contribution becomes 



VaU{x) = iV2/FTrdo(t){x) - ifliso[(T3, Uq] + 



(6) 



In the chiral Lagrangian Q), the mixed terms 
do4>{x)[a3, Uq] produce only boundary contributions and 
play no role here. Thus, to leading order in the e- 
expansion the fluctuation field ^(x) gives rise only to the 
kinetic energy term 



1 



where ^A — diag(m„,77id) is the quark mass matrix and 



d X -Tt d^_,4>{x)d^,4>{x), 



which decouples as in the theory with /iiso — 0. 

Collecting the remaining terms we see that the leading 
contribution to the partition function in the e-regime is 
the zero-dimensional integral 



Z{M;i^J,i, 



^jj ^iVF^,^l^TrlU,a3]lU-'^ ,a3] + ^SVTi {M^U+MU'') 



SU{2) 



(7) 



where we have dropped the 0-suffix on the group element U G SU(2). Projection onto fixed gauge field topology i^ 
[l5| is done by a Fourier transform, and amounts to the simple modification 



Ju{2) 



(8) 



r 



One sees that the leading-order contribution to the e- 
regime depends only on the scaling variables rhi = miYiV 
(where m.^ are the quark masses) and /i?^ = fif^^F^V. 
Both of these scaling variables are of order 1 in the e- 
counting. 

Effective partition functions related to (jHJ, and to its 
generalizations to more quark flavors of both kinds of 
statistics, have been studied in great detail recently [J, 
m m, III m 113 , and much has been learned about 
them. A particularly important feature for what follows 
is that such partition functions satisfy a series of exact 
relationships relating theories with different numbers of 
quark flavors to each other. The origin of the formalism 
lies in the theory of certain integrable systems, but we 
need here only the identities themselves, which are known 



under the names of Painleve and Toda lattice equations 
|28L |29| . These equations will be used to provide a non- 
perturbative definition of a replica limit, which in turn is 
needed to compute spectral correlation functions of the 
Dirac operator eigenvalues. 

We have organized this paper as follows. In Sec. II we 
reconsider the case of quenched QCD, for which a com- 
parison with lattice gauge theory simulations has already 
been presented [2|- In that paper the analytical results 
were stated without proof; here we provide the details. 
The main idea is to focus on a mixed spectral correlation 
function which is extremely sensitive to imaginary isospin 
chemical potential. In Sec. Ill we turn to the physically 
interesting case of two light quark flavors. The same two- 
point spectral correlation function is far more difficult to 



determine analytically. In a previous paper Q we briefly 
reported the final results, and showed how well they com- 
pare with lattice gauge theory simulations. The bulk of 
this paper, including all of Sec. Ill, is dedicated to the 
detailed derivation of just those results. Finally, Sec. IV 
contains our conclusions and an outlook on future work. 



II. QUENCHED THEORY 

In order to consider the quenched analogue of the sit- 
uation outlined in the Introduction we need to define the 
quenched limit on the effective field theory side. This 
issue was first resolved in Ref. J33 by means of a chi- 
ral Lagrangian living on a graded ( "supersymmetric" ) 
coset of spontaneous chiral symmetry breaking. An al- 
ternative, closer to the approach we shall pursue in this 
paper, relies on the replica method '31'|. We stress al- 
ready here that "quenching," be it by means of replicas 
or quark partners of bosonic statistics, is required even 
in the case of dynamical quarks if one wishes to com- 
pute spectral correlation functions of the Dirac operator. 



J 



pix+,X-;^^Mso)^{Y.s(x+-x^-AJ26(x--x^:^ 



Indeed, these methods are the only known approaches 
that allow access to the low-energy Dirac spectrum from 
effective field theory. 

The result of the quenched calculation was briefly 
stated in Ref. [J, which otherwise focused on the high 
numerical precision that can be reached for F^^ with the 
proposed method for measuring it. As a warm-up exer- 
cise for the Nf — 2 calculation we will here give the main 
ingredients behind this quenched result. We stress that 
the steps we follow are the same for both the quenched 
and dynamical cases; the only difference is that each step 
is simpler in the quenched case. We first define a two- 
point correlation function which is very sensitive to F^^. 
This correlation function can be obtained from a sus- 
ceptibility that we define and calculate in the effective 
theory in the e-regime. This last calculation, performed 
here using the replica method, is the most difficult part. 
It requires the use of generating functions that are ex- 
plicitly derived. 

The method of Ref. [21 is to consider the "mixed" two- 
point spectral correlation function of the Dirac operators 
D± defined in Eqs. ijl} and ©, 





En^+-^+ En^--^- 



.(m)^ 



(9) 



where the averages are performed over the pure Yang-Mills partition function. In order to reach the e-regime, this 
correlator is considered in the microscopic limit 






A generating function for the spectral correlation function jnj is the mixed scalar susceptibility, 



N / 1 1 

X(m+,m_; z^iso) = ( tr ^; tr ■ 



13+ + 771+ Z3_ + m_ 



tr- 



1 



13+ + mj 



tr- 



1 



D_ + 777_ 



(10) 



(11) 



As above, the averages are performed over the pure Yang-Mills partition function. (Note that this partition function 
is independent of to+ and 777_.) Written in an eigenvalue representation, Eq. (|ll|l becomes 



X(777+,777_;7Aiiso)=(5]-7;^^ E ^7;^ 



n ^A^ + ?7Z_|_ „^ %A_ + Vl\- 




(12) 



If one knows this function analytically, the spectral two-point function lO can be computed from the discontinuity 
across the imaginary axis (l^ |23| , 



/9(A+,A_;7/^iso) = ^Disc[x(777+,777_;^iso)]m+=iA+,m_=iA_ 
1 



47r^ e^0+ 

which is the inverse of the relation 

X(777+,777„;i/liso) 



73J lim [x(iA+ + e, iX- + e; ^iso) - x(«A+ - e, iA_ + e; ^iso) 



X(JA+ + e, iX- ~ e; ^iso) + x(«A+ - e, iX- - e; Aiiso)], 
p(A+, A_;7^iso) 



dXx- dX- 



{iX+ + m^){iX- + m_) 



(13) 



(14) 



In the replica formalism the mixed scalar susceptibility Hll() can be defined as [SJ 

x(TO+,m_;i/iiso) = lim -^d„i,dm_ logZ2„(m+, m_; i^iso), (15) 

where Z2j^{m^, ?n_; i/iiso) is the effective partition function of 2n replicated quark flavors. Half of these have degenerate 
masses m^ and chemical potential ifj-iso, while the remaining n flavors have degenerate masses m_ and chemical 
potential —i^iso- At the level of the fundamental theory, 

Z2"„(m+,m_;^^i,o) = AdA], det(i^+ + m+)" det(i^_ + m_)"e-^™(^). (16) 

The leading contribution in the e-regime is analogous to the SU(2) case discussed in the introduction, and reads J25l | 

Z2-„(m+,m_;z/iiso)= / df/ det(C/)'^e^'^?==Tr [r/,s][c/t,i3]+iTr (Mtcz+MC/t)^ ^^^^ 

where B — diag(l„, — !«) and M ~ diag(TO+, . . . , m+, to_, . . . , m_), in terms of the scaling variables 

m± = m±Y,V and (liso = ^iisoF^rVv. (18) 

We also use Z to denote the partition function in the e-regime, but to distinguish it from the general partition function 
(|16|l . we write it explicitly as a function of the scaled variables. Because these partition functions p7ll are linked to 
an integrable structure (of an otherwise unrelated Hamiltonian system), it turns out that all can be derived from just 
the n = 1 case |33 |. 

Jo 
by means of |25| 

22nim+,m_;ifiiso) = (^ ^ L(«-i) ^^^ [("i+5A+)''(m-9A_)'22''(m+, to_; i/ii^o)] ^_;^o_j_^^^_„_^ . (20) 

Here D„ is a normalization factor whose exact value need not concern us here. It is chosen so that the effective 
partition functions (|17|1 satisfy the Toda lattice equation (£)„ determines the coefficient of proportionality) (23 

~ o, - o , ^1/^- - •- \ A 2/- - \2^2n+2("^+'™-;*Aiso)2^2n-2(™+'™-;*Aiso) ,„i ^ 

m+dfn+m^dyyi_ logZ2„(m+, m„; i^iso) = 4n (m+TO_) , : — . (21) 

[22n("^+."l_;iAtiso)]^ 

Using this exact equation to define the replica limit, we obtain 

x(m+,m_;i/iiso) = 4m+m_Z2(m+,m_;i/iiso)-Z^2("T-+,"i_;i/iiso) ■ (22) 

Apart from the U(2) partition function IjHJl, now given by (|19f) . we also need 2^2('^+j^-;*Aiso)- This is the 
partition function of —2 fermionic quarks, i.e., 2 quarks of bosonic statistics. Like its fermionic analogue, the effective 
bosonic partition function is determined by the symmetries of the underlying QCD Lagrangian, now with bosonic 
quarks. In addition, in the bosonic case one must carefully enforce the convergence of the partition function |ig.l33l|. 
a problem that is absent in the fermionic case due to the nature of Grassmann integration. As the purely imaginary 
chemical potential does not change the hermiticity properties of the Dirac operator, this does not lead to additional 
constraints as it does for real chemical potential 2^, 27J . The bosonic partition function can be written as the integral 

Z'l2(^+,^-;Wso) = /^^4^det(Q)''e-i^?=°*'-[Q'^l[Q-\B]-itr(MtQ+Q-i^)^ ^23) 

J det Q 



where dQ 9{Q) /det Q is the integration measure on positive definite Hermitian matrices. Using the parametrization 

e*" cosh s e** sinh s 
e~*^sinhs e~''coshs 



„ */ e'"coshs e sinhs i ,„., 



where 

re [—00,00], se[— 00,00], te[— 00,00], 6'e[0,7r], (25) 



we find 



Z^2("i+i^-;*Aiso) =6^'^'=° / dtte^'^^^-°^^ K^{trh+)K^{tm-). 



At /tiso = this integral can be done analytically, 



Z'^2{m+,m-;ifliso = 0) 



(26) 



(27) 



and the result agrees with the expression at /Xjso = that was derived in Refs. 34, 35]. The main difference between 
this and the corresponding ferniionic result is the replacement of modified Bcsscl functions In{x) by Kn(x). This 
can be traced back to the non-compact integration range in (|23|l . which in turn follows from the symmetries and 
convergence requirements of the theory with bosonic quarks. 

With the above ingredients we immediately find the mixed scalar susceptibility from Eq. (|22|l . 



xim+,m^;ifliso) = Am+rh^ 



dt te^'^''-*'' I^ {tm+ )I^ {trh^ ) 



dtte'^i^^-"*^ Kt,{tm+)Kt,{trh^) 



(28) 



Taking the discontinuity as dictated by Eq. (|13|l . we finally obtain the desired quenched spectral correlation function, 






/ dtte^f''^^'"jM+)Mt^-) / dite-2AL*V,(te+)J.(ie-) 



(29) 



d<te2^-*V,(i^+)J,(te-) 



wM-^)<m:)-I>"'-''''"''^''*>'' 



(tC-) 



In the last line we have traded one non-compact integral for an integral over the compact interval [0,1]. This is 
convenient if one wishes to evaluate the expression numerically. 

Again it is useful to check the limiting case fj,iso = 0, where the above expression can be simplified. Using the 
orthogonality properties of modified Bessel functions on the interval [0, cx)], we have 



/ dttJ^{t^+)J^{t^-) = / dttJ^{t^+)J^{t^^)- / dttJ^{t^+)J^{t^+) 
Ji Jo Jo 

1 .., , . c+j.+i{^+)JA^-) - g^ j.+i(e-) J.(e+) 



(30) 



This allows us to rewrite the above expression for fiiso = as 



p,(^+,c_;z/ii,o = o) = 5(e+-e-)Y[j'(e+)-j.+i(^+)j.-i(^+)] 



ff2 _ ^2 \2 l^+'^'^+'i-iMJ^i^-) - S,-J,y+li(-)J,^i(^ 



(31) 



r 



which agrees with the known result [la. ITtI |2(]| | . When 
the two eigenvalues ^± coincide there is an explicit S- 
function contribution whose coefficient is given by the 
spectral one-point function. This is due to eigenvalue 
auto-correlation. The way that the finite imaginary 
isospin chemical potential /iiso "resolves" this (5-function 
contribution (because the eigenvalues ^_ and ^+ are now 
associated with two different Dirac operators) is quite 
spectacular. We show in Fig. ^ the quenched spectral 
two-point function (|29() with one eigenvalue arbitrarily 
fixed at f_ = 4 as a function of the other eigenvalue 



£_+ for /"4„ ^ 0, 0.001, 0.01, and 0.05. The pronounced 
peak around ^+ = 4 is precisely the remnant of the 5- 
function at ^_ = C+- It has been shown in Ref. 2] how 
this spreading out of the (5- function provides a method for 
determining the pion decay constant F^^ in lattice gauge 
theory simulations. 







-0.2 



4 



FIG. 1: The quenched two point-correlation function (|22 
with one eigenvalue arbitrarily fixed at 5_ = 4 as a function 
of the other eigenvalue 5+. The full curves correspond to 
A?so = 0.001, 0.01 and 0.05. The curve with the narrowest 
peak corresponds to the smallest value of fnao- Also shown 
(dashed) is the curve for fiiso ~ where the peak becomes a 5- 
function at ^_ = ^+ = 4. The (5-function has been suppressed 
from this figure. 



III. TWO LIGHT FLAVORS 

We now turn to the physically more important prob- 
lem: QCD with two light flavors. As stated above, we 
follow the same steps as in the quenched case. We first 
define the two-point correlation function that is very sen- 
sitive to Ftt- This correlator can easily be obtained from 
a susceptibility that we calculate in the e-regime of the 
low-energy effective theory using the replica method. To 
perform this calculation, we use a Toda lattice equa- 
tion which requires the introduction of several generating 
functions with different numbers of fermionic and bosonic 
quarks. 

We thus consider the correlation function 






(32) 



EH 



(") 



j:s(x--x 



(m) 



A, -A 



between the eigenvalues A^ and A_ of the anti-hermitian 
operators D^ and £)_, defined in Eqs. 1^ and (0), re- 
spectively. The average in H32|l is taken over the QCD 
partition function with two light flavors with masses ruu 
and 771(1, viz. 



J[dA]^ ■■■ dct{D++777u)det{D.+77id)e-^™'''^'^ 
'' /[dA]^ det(D+ + m„) dct{D^ + 77id)e-'5™(A) 

As in the quenched case, the e-regime is reached by taking the microscopic limit of this correlator, 

where 7fiu,d = mu,d^V. 
The scalar susceptibility, 



x(to+, TO_, m„, md; ifUso) = ( tr 



1 



-tr- 



1 



Z5+ + 771+ £)_ + 771- 

which in terms of the eigenvalues A-|_ and A_ can be written as 

1 'sr^ 1 



- tr 



1 



D+ + 771 < 



tr- 



1 



£)_ + m_ 



x{777+,7n_,mu,7nd;iHiso) = ( X! 



E- 



„ ZA^ + 7TI+ „j l\_ + m^ 



E 







■ \ (™) 

IA_ + 771^ 



(33) 



(34) 



(35) 



allows us to calculate the correlation function from the discontinuity across the imaginary axis since the QCD partition 
function does not depend on r7T,±: 



p{X+,X-,mu,7nd;i^J,iso) = — ^Disc[x(m+,7TT._,77i„,7TT,<i; A^: 



iso J Im-u—iX-i- .m_ —iX- 



= -j-^ lira [x(«A+ + e, iA- + e, r77„, 771^; A^iso) - x(«A+ - e, i\- + e, 77i„, tti^; ^i^o) 



47r2 £-,0+ 



-x(»A+ -1- e, i\^ - e, r77„, tti^; /^iso) + x(*-^+ - e^ «^- - ^i w„, Tn^; /Xiso)]. (36) 



The inverse of this relation is 

A. The susceptibility from the rephca limit 

In order to obtain the susceptibiHty we again employ the replica method, writing 

x(TO+,m_,m„,TOd;iA'iso) = 1™ ^j^™, a,„_ log Z2„ 2 (™+i '7i_, m„,m<i; i/ijso)- (38) 

The generating functions Z2„ 2 have 2n replica flavors in addition to the two flavors of mass niu and md, 

Z2n,2i^^+^'^~iiT^u,md;ifiiso) = /[d^]i. [det(D+ + m+)det(I?_ +m_)]"' 

X det(i:'+ + rriu) det(i)_ + md)e~^™'-^'> . (39) 

Note that half of the replica flavors have mass rn+ and chemical potential i/iiso while the other half have mass rn_ 
and chemical potential —i^iso- 

In the e-regime the leading contributions to the partition functions -E2„ 2 again satisfy Toda lattice equations. To 
obtain the correct replica limit n ^ in Eq. H38() we make use of this Toda lattice equation [23 : 



■m+drh_^rh-dm_ \ogZ2^2{'i^i+,m-,mu,md]ifi.iso) (40) 

= 4n {rh+ifi-) 



2,- - ^2^2n+2,2(™+>™->'^n,md;^/iiso)2^2„„2,2('T^+,TO_,TO„,md;^Aiso) 



[-2^2''«,2("^+' "1-7 "i«> "^d; *Aiso)]^ 

Taking the replica limit we arrive at 

,. . . . .„ ^ .- - Z|'2(TO+,rn_,m„,mrf;i/iiso)2^2''-2('^«i"id|"i+,m_;i/iiso) ,^_,, 

X(m+,TO_,TO„,md;j^iso) = 4TO+m_ — ■ r— — -^ :-^. --^ , (41) 

[Z2(m„,md;i/iiso)J 

where 22„ 2 has been defined in Eq. (|39|) : Z2 {mu,rhd;i fiiso) — -2g 2(m+, m_, tti^jTO^; /iiso) is the partition function 
with zero replica flavors; and Zj _2(TOti,md|TO+, m_; i/tiso) = Z^2 2(''^+j '^-i™")"^^; Aiso)- 

The discontinuity of x across the imaginary to-|_ and m- axes gives the dynamic correlation function, 

,^ ^ . .... ^ ^ Zl2{'i^+,'^£.-,'n^u,md;iiJiso)Z^-2{mu,md\i^+,i£.-;ifiiso) ,,^, 

p.(^-.,?-,-.,-.;^Miso)=e+e— [Z2^(^„,^,;.^,„)]2 • (42) 

We therefore need to calculate three different generating functions: two with fermionic quarks only, and one with 
both fermionic and bosonic quarks. 

B. Computing the fermionic partition function 

In the eifective theory in the e-regime, the generating functions with n > involved in the calculation of the 
susceptibility are given by 

Z2^„2(m+,m-,m„,m,;zAiso) - / dU det([/)''ei'^ioT'-[^-^][^nBl+iTr(Mtf/+MC/t)^ (43) 

with 

B = diag(l„+i,-l„+i) and Af = diag(m+, ..., m+,TO„, ?fi_, ..., m_, m^). (44) 

The partition function Z2 (tok, m^; J/iiso) has been calculated in Eq. (|19|1 above. In addition, 

Z2 2{iT^+i'iTi-T^u}^dj'i'P-iso) can be obtained from the calculation presented in Ref. 24] by changing the sign 
of /i?,Q. We have thus 



Z2.2(m+, m_, TO„, fiid; ifiiso) 



(to^ — ml^){m^ — rri^) 



Z^(m+,TO_; j/liso) Z^{m+,md;ifiiso) 
Z^{mu,rh-;ip,iso) Z^{mu,md;ifj'iso) 



(45) 



C. Computing the supersymmetric partition functions 

The calculation of Z2 _2(»ti„, TOd|m+, m_; ifiiso) requires us to perform an exact integral over the supergroup Gl(2|2) 
with two fermionic and two bosonic quark flavors. This is a rather lengthy analytical calculation. We start from the 
effective generating function in the e-regime which is given by 



^Gl(2|2) 



r 



(46) 



where B = diag(l; 1, — 1, — 1), and the mass matrix is 
given by 



M = 



M+ 
M_ 



(47) 



with Af+ = diag(TOi,,m+) and Af_ — diag(md,m_) pjq . 
At /ijso = 0, this is the same as the generating func- 
tion used in ^^. Notice that in the complete effective 
partition function at ^iso = there is an invariant op- 
erator that contains the term i9^Str InU = Str If^d^U 
[isL [13 . This operator is obviously absent from the zero- 
momentum part of the partition function. At non-zero 
/iiso, this operator is also absent from the zero-momentum 
part of the partition function since [3g StrC/~^VoC/ — 
Str U-'^doU, and VqU = OqU - ijHsolB, U] . We can un- 
derstand this from physics as well: The isospin singlet 
does not couple to isospin chemical potential. 



1. Parameterization of the Goldstone supermanifold 

In order to calculate the exact supergroup integral 
(|46|l . we have to parameterize the Goldstone manifold. 
We use the same factorizing parameterization as in 
Ref. IB: 



U 



wi 

W2 



\/l — WW 

—w a/T 



w 



WW 



wi 

W2 



(48) 
where w;i,2 £ Gl(l|l) and w, w e Gl(l|l). This pa- 
rameterization leads to rather simple expressions for the 
integrand of (|1^ : 



1 



SdetU = Sdet(u;^ wj), 



StrMiU + U-^) = Si+S2, 



and 



with 



iALStr[i3,C/][B,C/-i]^5^, 



(49) 
(50) 

(51) 



iSi = -Str [Af+ (wiVl — WW wi 



+w^'VT 



WW w 



r^)], 



(52) 



52 = -Str [M_ (w2\/r 



WW W2 



+W2^VT 



WWW2^)] , 



and 



Sf, = -2/if,„Str [ww] . 



(53) 
(54) 



For w^ 2 G Gl(l|l), we use the same parameterization as 
the one used in Ref. .19], 



w,- 



^Vi, 



where 

A,: 







and Vi — exp 



ai 

A 



(55) 



(56) 



with ipi g (— tTjTt) and Si G (—00,00); Ui and Pi are 
Grassmann variables. The main advantage of this pa- 
rameterization is that its Berezinian is equal to 1, as was 
shown in Ref. [l9j. Thus we parameterize the matrices 
wi,2 € Gl(l|l) so that wf = h-iVi, i.e.. 



(1 + \c,a,Pi) e^^'/2 



(1 



aiUi 



,(57) 



where 



Ji'i 



e»«A>/2 -|_ eSi/2 



(58) 



1 

Ci = — 
2 

di — - 
2 



pji/;,/2 



./2)2 



aii'i 



piV./Z 



./2)2- 



For w^w £ Gl(l|l) we use the same polar decomposition 
as in |23: 



w = vSu 



and w = uSv 



(59) 



where S and S are 2x2 diagonal supermatrices with 
commuting elements given by 

S = diag(sin6'e*'',zsinh0e*''), (60) 

S = di'Ag{sm0e-'P,ismh(l)e-'''), (61) 



with 9 e (0,7r/2), cf) € (0,oo), and p,a € (— tt, tt), while 
u,ve U(1\1)/[U{1) X C/(l)] are given by [13 



u = exp 



c 



and V — exp 



e 



(62) 



X ; ""- - -^' V ^ 

With these parameterizations, we obtain 

SdetU = Y[e''^'-'\ (63) 

i 

In addition, the supertraces in Eqs. itCT) . (|55|l . and itS^ 
are given by 

Si = ^Str [wC^;~i (wiM+wi + w^Hl+w^^)] , (64) 



52 = -Str [liCw-i (w2A/-W2 + w^Hl^w^^)] , (65) 
and 



5^ = 2Af3oStrC^ 



(66) 



where C = VT^-SS ~ diag(cos6',cosh0). The advan- 
tage of our parameterization is that the integral over the 
supergroup Gl(2|2) explicitly contains two independent 
integrals over Gl(l|l), which are simpler to compute an- 
alytically. 

In order to perform the group integral we need to deter- 
mine the integration measure that corresponds to our pa- 
rameterization. The parameterization of the Goldstone 
manifold (|48|l is of the form 



U = WTW. 



(67) 



As was shown in Ref. [2j|, the measure factorizes into a 
product of one factor that depends only on W and one 
factor that depends only on T, 



dfi{U) = w^'^ dw 



1 "'2 



' dw^ T-^ dT 



For wf 2 & Gl(l|l), we have used the same parameteriza- 
tion as in Ref. ,li^] . In that paper, it was shown that the 
Berezinian of this change of variables is equal to 1 , and 
thus that 



w~ dwf = d^i dsi dai d(3i . 
We therefore find that 

fi{'Wi) dwi = dsi dipi dai d(3i 



(69) 



(70) 



Finally, the parameterization we use for T is exactly the 
same as the one used in Ref. [23 ■ The measure is given 
by 

T^ dT = ij,{w , w) dw dw 

— 2^-^d6 d(f>dpda dC,dxd^dri. 



(cos^ ~ cosh (f>)'^ 

2. Efetov-Wegner terms 



(71) 



= p,{wi)dwi p{w2)dw2 p,{w,w)dwdw. (68) 



As for any supersymmetric integral, extreme care has 
to be taken with the singularities that might be intro- 
duced through a specific parameterization of the inte- 
gration supermanifold and the corresponding measure. 
The singularities of the measure affect the supersymmet- 
ric integral through the so called Efetov-Wegner terms. 
(See for example Refs. [ll,|33| for a discussion of Efetov- 
Wegner terms.) The measure p.{wi) dwi does not con- 
tain any singularity in the variables Si and ^i, and 
there are no Efetov-Wegner terms related to our param- 
eterization of wi^2- On the other hand, the measure 
T~^ dT — p,{w,w)dwdw is singular when 9 = i(j). We 
therefore expect Efetov-Wegner terms in this case. The 
method used in [33| can be straightforwardly applied to 
compute the Efetov-Wegner terms related to our param- 
eterization of T. With our parameterization, including 
the Efetov-Wegner terms, we find that the generating 
function 2^2-2 i^ given by 



Z2 _2imu,rnd\m+,m^;ip,iso) 



lim / d9 d(f) dp da dC, dx d^ drj 



where 



(27r)3 e^Q 



(cos^ 9 — cosh 



^A\2 



(72) 



{vc - e) + S{vc - e)vn + -5'{vc - e)vl 



11X26^^^°^'^°^^ S-cosh^ (f>) 



Vc = sin^ 9 + sinh^ 

Vn = (sin^ 9 - sinh^ 0) {(x + C??) + 2i sin 9 sinh cj) (e'^P-"^ r/C 
-H2(sin^ 9 + sinh^ (p)Cx£,V 



_ g-«(p-'T) 



ex 



2(sin"'6i-|-sinh^(/))^CxC77 



(73) 



and 



1 pTT /'OO 

I. = — / d^, / ds, I da, dp, e^^+-("^^--^) 
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(74) 



The normalization is chosen so that Z2 _2(^u'™d|?^«,™d; ^Aiso) = 1- 



3. Analytical result for the supersymmetric partition function 

We are now in position to expUcitly compute the partition function (|46|l . We first compute Ti given by Eq. 174|l . 
The integral over a^ and Pi is readily obtained by expanding exp5i to first order in a; and Pi. We get 



Xi = J'(m„,m_(., cos 6*, cosh 0) + ^r7/C(TO„,m+, cos 6*, cosh (/)), 
I2 — J^ {rh d, rfi -, cos 0, cosh. <p) + (^xl^i'>TT'd j'fiT'-j'^os 9, cosh (p), 



(75) 



with 



1 r pi'(JV'-s) 

J(x,y,t,p) = / dVdse*^™''^"''^™'^' ^ 

^^ '^' '^^ 16^7 ^ 1 + cosh^ 



2tx cos ^/i + 2py cosh s + py cosh 



3s — itp 



{2px — tx ~ py + 2ty) cosh 



s -\- iip 



(76) 



'tx cosh 



s — Siip 



and 



lC{x,y,t,p) 



t_p 
647r 



dipds e 



tx cos ^—py cosh s iy{i^ — s) 



X [px'^ — ty^ + x{2 cos ^ + tx cos 2^) + y (2 cosh s — py cosh 2s) 
+2i(t — p)a;?;sin'0 sinhs] . 
Hence the partition function <|72(l can be written as 



(77) 



Z2_2{mu,md\m+,m_-,ifJ.iso) = 1™ dt f dp 



tp 



e^O /n 1 1 [t'^-p^Y 



^d{p^ ~e~e) 



xJC{rhu,m+,t,p)}C{rhd,rh^,t,p) 



"'0 



+ lim I du I dv—^ ^-^7 



eA.iu-+-')S{u^ + v^-e) 



(78) 



X Hu^ — w^) /C f m„,m+, vl - u^, v 1 + u^ j J' (TOd,m_, v 1 - u^, -y/l + ti2 j 
+J f m„,m+, V 1 — u2, ^/l + t;2 j /C (TOd,TO_, v 1 — u^, V 1 + I'^j 
+ (m^ + v^)J (rhu, m+, v^T^^, V^ + ^A J (™d, m-, Vl - u^, ^1 + i'^) } , 
where t = cos 9, p = coshc/), u = sm9, and v — sinh0. We can rewrite the second integral as 

|.7r/2 i-l/sinq _^ ( 

lim / dq dr sm2qe~^^^=°'^S{r — e)< {cos^ q ~ siii^ q) 

^-oJo Jo [ 

X /C I m„,TO+, \/l + rcos2(j, y 1 ^ rsin^q I >/ ( rhd,rh-, \/\ + rcos'^q, y 



1 — r sin g 



+1/ I w„, m_|_, v 1 + r cos2 g, \/l — rsin 5 I A^ I rhd,rh-, v 1 + rcos2 g, \/ 



' 1 — r sin g 



+J ( mu,m+, v^l + rcos2g, y 1 - r sin^ q 1 J ( rridjm-, \/l + rcos^q, y 
J{mu,m+, 1, l)J{rhd,rh^, 1, 1), 



1 — r sin g 



(79) 
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where u = ^/r cos q and v — \fr sin q. This gives 



Z^_^2{™'u, TOd|m+, m^;ifliso) 



tp 



o'V, ''\t^-P^y 



,2AL(*'-p')/ 



/C(m„,m+,t,p)/C(m.>.,t,p) 
+J^(m„,m+,l,l)J^(md,m_,l,l). (80) 



Finally, we have to calculate ^{x, y, 1, 1) and IC{x, y, t,p). The result is 

1 



J{x, y, 1, 1) - 2 ^^ [^-i(^) + ^-+i(^)] ^-(y) + yl-(^) [K.-iiy) + K,+i{y)]} ., 



and 



IC{x,y,t,p) = {x^ - y^){t^ - p^)I,{tx)K,{py). 
The partition function can finally be written as 
Z2^_2{mu,fhd\m+,m^;ifliso) = [(m„)^ - (?7i+)^] [(?7irf)^ - (m_)^] 

X / dtte2^-*/^(im„)/^(imd) / dppe'^^^-"^ K^{p'm+)K„{prh^) 
1 



(81) 
(82) 



(83) 



+ - {rhu [h-iirhu] + /i.+i(m„)] K^{m+) + m+I^{rhu) [K„^i{m+) + ^4:^+1(771+)]} 
This can be written in a more compact notation as [cf. Eqs. H19|l and (|26|l ] 



Z2_2{mu,md\m+,m-;illiso) = 



{ttiI -m\)Z!^{mu,rhd]ifi'\so) Z'(_^{md\rh^) 

{ml - m?. )Z^2 (^+ : "1- ; ^Aiso) 



(84) 



-ZJ'__i(m„|m+) 
where the /ijso-independent graded partition function, 

-2r,-i("iu|TO+) = - {rhu [h-i{mu) + /^+i(m„)] K^{m+) + m+I^{mu) [K^^i{m+) + K^+i{m+)]} , (85) 

was calculated in Ref. 19|. 

Note that Zj _2('7i+,?7i_|m4-,TO_;i/iiso) = 1 as it should. Furthermore, upon expanding 

-E2 _2(w„,md|TO+, m_; i/tiso) to leading order in rriu — m+ and rud — rii-, the quenched correlation function 
(|29ll is recovered using the supersymmetric method. 



D. Final result 



We can thus finally compute the two-point correlation function H42|l . using the analytically calculated generating 
hmctions ^^, (gHJ, and ^^. The result is 



1 ^ ^ -,-2 

dtte^^^^-^^h{tmu)lAtmd) 



Ps{£.+ ,^-,mu,rhd;ifiiso) = S+C- 

x( I dtt e2'^?-*Vj.(i ^+)JAt C-) / dt t e2^-°*'/^(i m,J/^(t m^) 
dt t e2^-*'/,(t 7fiu)JAt C-) / dt t e^^-*' J,(i ^+)L{t rhdyj 



X / dtte^^'^'''''l^{tmu)I^{tmd) dt i g-^^'-*' J^(U+)^<.(i ^-) 



(mu^+i("^u)J.y(g+) +C+'^.^+i(C+)^t/("^u))("^<i-?"iy+i(?Tid)J.y($-) +g-J^+i(g-)/t/(mrf)) 

I 



(86) 



For a numerical evaluation it is advantageous to rewrite in Eq. (|29|l . 
the non-compact integral appearing in the fourth line as 
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We have performed various checks on this result. For 
example, at /tiso = it correctly reduces to the two- 
point microscopic correlation functions at zero chemical 
potential ^^ [13 ■ We have also verified that it reduces 
to the quenched result H29() in the limit where both to„ 
and rud are sent to infinity, as is required by decoupling. 

As in the quenched case the correlation function at 
Miso — has a (5-function at equal arguments. When 
/iiso is nonzero this i5-function becomes a peak in the 
correlation function around equal arguments, as shown 
in Fig. 13 



0.2 



0.15 - 



\i=0.0 

fi,,;F;V=o.o5 




-0.05 - 



-0.1 



^. 



FIG. 2: The two-point correlation function I86L with one 
eigenvalue fixed at ^^ = 4, as a function of the other eigen- 
value ^+ . The masses of the two dynamical flavors are chosen 
to be degenerate at the value m„ = rhd — 5. The dashed curve 
shows the correlation function for fiiao = and the solid curve 
corresponds to fJ.fsoF^'^ = 0.05. The (5-function that appears 
at ^+ — (,- — 4 when /iiso = has been suppressed from this 
figure. 



the physical case with two dynamical light quark flavors. 
In two previous articles, it has been demonstrated that 
these formulas for the correlation functions lead to an 
efficient way to determine i^^ on the lattice 0, y| • 

Our calculation made extensive use of the replica 
method. In this approach, the correct answer is reached 
via Toda lattice equations, which relate theories with dif- 
fering numbers of quark species to each other. Taking the 
replica limit required the computation of several partition 
functions with varying numbers of bosonic and fermionic 
quarks, which led us to compute some non-trivial parti- 
tion functions. Our strategy thus required the calcula- 
tion of exact group integrals over both graded and non- 
graded Goldstonc manifolds. It must be noted that our 
calculation would be technically very tedious if it were 
carried out exclusively by the so-called supersymmetric 
approach: It would require an integration over Gl(4|2), 
a highly complicated task. The advantage of the Toda 
lattice equation is that it reduces the complexity of the 
calculation by spreading the difficulty over several parti- 
tion functions. 

As a tool to extract F^^ from dynamical lattice sim- 
ulations it is of obvious interest to derive the correla- 
tion function in a partially quenched theory where the 
chemical potential is set to zero for the physical dynam- 
ical quarks. This would allow the use of existing gauge 
field ensembles, generated at zero chemical potential, in 
a determination of F^^ . The analytical expression for this 
partially quenched correlation function has proved to be 
challenging. We hope that the detailed calculation pre- 
sented here may be of help in the future calculation of 
such quantities. 
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